An Exact Solution for Spin and Charge Correlations in Quantum Dots: 
The Effect of Level Fluctuations and Zeeman Splitting 



I.S. Burmistrov 1 , Yuval Gcfcn 2 and M.N. Kiselev 3 
1 L.D. Landau Institute for Theoretical Physics RAS, Kosygina street 2, 119334 Moscow, Russia 
Department of Condensed Matter Physics, The Weizmann Institute of Science, Rehovot 76100, Israel and 
3 International Center for Theoretical Physics, Strada Costiera 11, 34014 Trieste, Italy 

(Dated: January 24, 2012) 

The inclusion of charging and spin-exchange interactions within the Universal Hamiltonian de- 
scription of quantum dots is challenging as it leads to a non-Abelian action. Here we present an 
exact analytical solution of the probem, in particular, in the vicinity of the Stoner instabilty. We 
calculate the tunneling density of states and the spin susceptibility. We demonstrate that near the 
Stoner instability the spin susceptibility follows a Curie law with an effective spin. The latter de- 
pends logarithmically on temperature due to the statistical fluctuations of the single-particle levels. 
Near the Stoner instability the tunneling density of states exhibits a non-monotonous behavior as 
function of the tunneling energy, even at temperatures higher than the exchange energy. This is 
due to ehnanced spin correlations. Our results could be tested in quantum dots made of nearly 
ferromagnetic materials. 

PACS numbers: 73.23.Hk, 75.75.-c, 73.63. Kv 



I. INTRODUCTION 

The study of quantum dots (QDs) is at the cutting 
edge of modern condensed matter physics. The introduc- 
tion^ of the universal Hamiltonian (UH) made it possi- 
ble under not severe assumptions to describe a variety of 
QDs by means of an effective zero-dimensional Hamilto- 
nian with a few physical parameters. In particular, the 
UH allows one to simplify the intricate electron-electron 
interactions within a QD in a controlled way. 

Within the framework of the UH electron-electron in- 
teraction is represented as the sum of three spatially inde- 
pendent terms: charging (singlet particle- hole channel), 
spin-exchange (triplet particle- hole channel), and interac- 
tion in the Cooper channel. The latter is responsible for 
superconducting correlations in QDs. In what follows, we 
shall assume that the Cooper channel is suppressed, e.g., 
by the orbital effect of a weak magnetic field. The charg- 
ing term is responsible for the well-known phenomenon 
of Coulomb blockaded It is broadly known that in the 
presence of significant ferromagnetic spin-exchange in- 
teraction, bulk systems can undergo a Stoner transition 
from para- to ferromagnetic materials. In the case of 
QDs physics is richer^ One distinguishes three regimes 
of behavior as function of the increased strength of the 
ferromagnetic exchange interaction (J > 0): i) param- 
agnetic (the total spin in the ground state is zero); ii) 
mesoscopic Stoner regime (finite total spin in the ground 
state whose value increases stepwise with the exchange) ; 
and iii) thermodynamic ferromagnetic phase (the total 
spin in the ground state is proportional to the volume of 
a QD). The mesoscopic Stoner regime disappears in the 
thermodynamic limit: 8 — > 0, where 6 is the mean spac- 
ing between single-particle energy levels in a QD. The 
mesoscopic Stoner regime is sensitive to the statistical 
fluctuations of single particle levels^ and to the presence 
of the Zeeman splitting^ The former enhances the total 



spin in the ground state whereas the latter suppresses the 
mesoscopic Stoner instability. To take both effects into 
account simultaneously calls for a full-fledged quantum 
mechanical treatment of the problem. 

At first glance, the UH with charging and spin- 
exchange interaction terms is easy to solve. All its 
three parts (free electron term, charging term, and spin- 
exchange term) commute with each other. It allows one 
to work in a basis of states classified by the total number 
of electrons and the total spin^r— However, this approach 
requires calculation of Clebsch-Gordan coefficients which 
is not an easy task. In this way Alhassid and Rupp* have 
found an exact solution for the partition function in the 
absence of Zeeman splitting. Elements of their analysis 
were then incorporated into a master equation analysis of 
electric^ and thermal^ transport through a QD at low 
temperatures T < 6. Independently, a study of electron 
transport through a QD for the same temperature range, 
T < 5, was made by Usaj and Bar anger & Their analysis, 
accounting for the charging and exchange interactions, 
was based on a master equation approach as well. 

More traditional way to tackle an interacting prob- 
lem is to employ the Hubbard-Stratonovich transfor- 
mation. The latter reduces the interacting problem to 
the problem of non-interacting electrons in the presence 
of dynamical (time-dependent) external field (see e.g., 
Ref. [l(| ■ For the case of the UH with the charging inter- 
action only such dynamical external field can be removed 
by suitable gauge transformation of fcrmionic opera- 
tors. This method was first employed by Kamcnev and 
Gefenii and resulted in the exact solution of the prob- 
lemi 12 i 13 The gauge transformation can also be treated by 
means of stochastic (Langevin) equations. Correspond- 
ing stochastic bosonization approach based on construc- 
tion and solution of Fokker-Plank equations was used in 
Ref. [3| to solve exactly the UH with the charging and 
Ising-spin-cxchangc interactions. 
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In the case of isotropic spin-exchange interaction the 
Hubbard-Stratonovich transformation results in the ef- 
fective action which describes electrons with a dynam- 
ical time-dependent magnetic field 6(t) acting on their 
spinsJ^ Solving this non-Abelian effective action is an 
intricate problem. One needs to tackle time ordered ex- 
ponents of the form 

Tcxp dt' d{t')s{t') \ , (1) 

where s represents an electron spin and T is a time or- 
dering operation. To avoid the problem of the time or- 
dering the UH with anisotropic spin-exchange interaction 
was considered, and perturbation expansion around Ising 
point was performed.— 

Here we present an exact analytic algorithm to tackle 
the challenging problem of the UH with charging and 
isotropic spin-exchange interactions in the presence of 
Zeeman spliting. To solve the problem of a time or- 
dered exponent we employ here a Wei-Norman-Kolokolov 
(WNK) transformation ! 16 i 17 Wei and Norman^ address- 
ing the problem of a quantum spin subject to a prescribed 
classical time-dependent magnetic field, have elegantly 
shown that by preforming a non-linear transformation 
from 6 to a set of new variables, the time ordered expo- 
nent (HJ) can be written as a product of three ordinary 
exponents (cf. Eq. (|24p). Even so, that problem could 
not be solved in a closed form in general. To find rela- 
tion between new variables and one has to solve Riccati 
equation. Although the problem of a dynamical magnetic 
field seems to be even more intricate, in fact, as it was 
shown by KolokolovjiZ it is simpler. In this case, one 
needs to know the Jacobian for the non-linear transfor- 
mation only. Then, the functional integration over new 
variables can be performed exactly. 

We thus present here exact analytic resuts for the par- 
tition function (cf. Eq. ([57]) ) and the tunneling density 
of states (cf. Eq. (|76p) for the UH with charging and 
isotropic spin exchange interactions in the presence of 
Zeeman spliting. We emphasize that our results are valid 
for arbitrary parameters of the UH. 

In the mesoscopic Stoncr regime, near the Stoner insta- 
bility, 5 — J <C 6, our general results can be drastically 
simplified. We find that in a wide temperature range 
J" « T < 5 J I (5 — J) the average zero-field spin suscep- 
tibility behaves according to the Curie law with a large 
effective spin which depends on temperature logarithmi- 
cally (cf. Eq. ([TUT]) . The latter is the effect of statisti- 
cal fluctuati ons of single- particle levels. A tiny magnetic 
field B ~ y/ JT{1 — J/6)/g/j,B is enough for the average 
spin susceptibility to become temperature independent 
Fermi- liquid like (cf. Eq. (|111[) ). Here g and [is stand 
for the g-factor and the Bohr magneton, respectively. We 
find that enhanced spin correlations resulting in a large 
total spin in the ground state of a QD in the mesoscopic 
Stoner regime near the Stoner instability, 5 — J <C 6, can 
be observed as additional (to Coulomb blockade) non- 
monotonic behavior in the tunneling density of states 



(TDOS) at high temperatures S < T < 6 J/ (6 - J). 
Magnetic field suppresses the spin-related non-monotonic 
behavior of the TDOS. We mention that some of the re- 
sults were published in a brief form in Ref. [l8|]. Our 
main new results concern the effect of Zeeman splitting 
(cf. Eqs. (J7JJ, (HSU) ) and of disorder (cf. Eqs. ®, (EH) , 

(HQ), (UH)). 

The physics discussed in current work can be tested 
in QDs made of materials close to the thermodynamic 
Stoner instability, e.g., Co impurities in a Pd or Pt host, 
Fe or Mn dissolved in various transition-metal alloys, Ni 
impurities in a Pd host, and Co in Fe grains, as well as 
new nearly ferromagnetic rare-earth materials^ - — Pos- 
sibly, the intriguing magnetic behavior observed recently 
in Pd nanoparticles capped with different protective sys- 
tems^ is related to the physics of mesoscopic Stoncr 
regime. 

The outline of the paper is as follows. In Sec.|H]we in- 
troduce the UH, subsequent imaginary time action, and 
partially disentangle charge and spin degrees of freedom 
in the problem. In Sec. lIIII we introduce the WNK trans- 
formation to solve the problem of spin dynamics. In Sec- 
tions lIVI andfVlwith the help of the WNK transformation 
we derive exact analytic expressions for the grand canon- 
ical partition function and for the TDOS corresponding 
to the UH in the presence of Zeeman splitting and for 
a given realization of single-particle levels. The analysis 
presented in Section |VT|) incorporates the effect of disor- 
der. In Sections [VI AHVI Cl we present a rigorous analysis 
of the effect of level fluctuations on the spin susceptibil- 
ity. The latter is modified by disorder in a strong and 
significant manner. We refer the less initiate reader to a 
semi-qualitative derivation of our results for the average 
spin susceptibility (Sec. IVII"[) . In Sec. I Villi we discuss the 
dependence of the TDOS on energy, temperature and 
magnetic field. The effect of level fluctuations on the 
TDOS is semi-qualitatively discussed in Sec. HXl We con- 
clude the paper with summary of the main results, and 
brief coments of the amenability of our predictions to 
experimental tests (Sec. |X[) . 



II. FORMALISM 

A. Universal Hamiltonian 

We consider a quantum dot of linear size L in the so- 
called metallic regime, whose dimensionless conductance 
.gxh = Exh/S 3> 1. Here i?Th is the Thouless energy. We 
account for the following terms of the universal Hamilto- 
nian 2 

H = Hq + He + Hs, H = ^2 (2) 

a,tr 

Here, e ajCT denotes the spin (a) dependent single particle 
levels. In what follows, we shall assume that the magnetic 
field B is applied and t a . a = e a + g^sBa jl. 
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The charging interaction 

H c = E c (n- Nof 
accounts for the Coulomb blockade. Here 



(3) 
(4) 



is the particle number operator, and Nq represents the 
background charge. The term 



He = -JS 2 



(5) 



represents spin interactions within the dot. Here 

S = Sa = 2 51 a \^ a <y<y' a a,a' (6) 

a acrcr' 

denotes the operator of the total spin of electrons on the 
dot, with the components of er comprising of the Pauli 
matrices. 

We stress that we do not consider the interaction in 
the Cooper channel in the UH (J2J). For QDs fabricated 
in 2D electron gas the interaction in the Cooper chan- 
nel is typically repulsive and, therefore, renormalizes to 
zeroi In the absence of spin-orbit interaction the paral- 
lel magnetic field does not affect the orbital motion of 
electrons in a QD (we neglect the effect due to a finite 
width of 2D electron gas). In this case the statistics of 
single-particle energies e a can be described either by the 
orthogonal Wigncr-Dyson ensemble (class AI) or by the 
unitary Wigncr-Dyson ensemble (class A) i 23 ' 24 The lat- 
ter is achieved in a weak perpendicular magnetic field 
B± > B c = $ /(L 2 ,/<7Th) where <J> denotes the flux 
quantum. In the case of 3D quantum dots realized as 
small metallic grains, the interaction in the Cooper chan- 
nel can be attractive, giving rise to superconducting cor- 
relations. In this case we assume that there is a weak 
magnetic field B > B c which suppresses the Cooper 
channel. Therefore, the level statistics is described by 
the unitary Wigner-Dyson ensemble. 

The imaginary time action for the system reads 



J Cdr = J [^* Q (d r + I i)* a -H 



dr, (7) 



where [i is the chemical potential, j3 = 1/T, and 
we have introduced the Grassmann variables ty a — 
("tpaf: ipai) T , = (V'afj V'aj,) to represent electrons on 
the dot. 

Employing the Hubbard-Stratonovich transformation 
leads to a bosonized form 



g[i B B a 3> 
o T - e a — a z + /.i + t<p-\ — 



* 2 <f> 2 

+ — iN 6 

4J 4E C ° V 



(8) 



where <f> and <fr are scalar and vector bosonic fields respec- 
tively. The SU{2) non-Abclian character of the action 



poses a serious difficulty. In the presence of the charging 
interaction only (Abelian U(\) case) the problem can be 
solved by performing a gauge transformation* 1 ^— For 
the case of the charging interaction and spin-exchange 
interaction of Ising type (Abelian t/(l) x U{1) case) 
the problem also can be solved by a gauge transforma- 
tion] 14 ' 15 In the non- Abelian U{\) x SU(2) case, we start 
from performing a gauge transformation in the charging 
sector only. 



B. Partial disentanglement of spin and charge 

Our aim is to compute the grand partition function 
Z = Trexp(— 0H + /i/3n) and Green's function in the 
Matsubara time domain 



Tra Q , jCri (Ti)aJ (i0 . 2 (T2)e 



G a ^i,a 2 (TIj t 2) — —T T 

(9) 

Here, T T denotes a time ordering operation along Mat- 
subara time. In the Lagrangian formalism, the Green's 
function can be written as 

G (n,r 2 ) = ~% f P[*,*, ( /»,*]* Q (r 1 )^ C[ (T 2 ) e - s -, 



z ■= / X>[*,*,0, *] e - 5tot . 

Let us split the field </>(r) as 



(10) 



0(t) = 4>{t) +2™T + 4 / dT(j)(T)=Q, (If) 

Jo 

where the static component of 4>(t) obeys inequality 
1 00 1 < ttT. Then the part Mr) + 2-KmT of 0(r) can be 
gauged away (see Refs. pLlI — il5j ] for details). The Green's 
function becomes 

ttT 

G a (r 12 )= [ |%^D(r 12 ,0 o )^(ri 2 ,0o), (12) 



2ttT Z 



-TTT 



Z = 



ttT 



ttT 



2ttT 



(13) 



where t\i = t\ — t 2 . The so-called Coulomb-boson prop- 
agator reads 

D(t, 4> ) = e -Bc|r|(x-|r|/« e *(^+-) 

x e -l3E c (k-N 0+T /f3)\ / 14 x 

The Green's function <? a (ri 2 , 0o) on the right hand side 
of Eq. (|T2)| is defined as 

G a {.Ti2,4>o) = -T T / ^JT\ *a( T l)*Q(T 2 )e , 



(15) 
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Here the action 

rP 



S 



/ dr^c 
Jo 



d T - £q + p + ifio + 



4J 



(16) 



It can be formally rewritten as 



S = 



-H 



dr, H=Ho + H s , (17) 



where Ho is given by Hq (Eq. ([2])) in which e ai<T is re- 
placed by e a>CT = e Q!(T - /i + i<fo: 



'Ho — /J £q 



(18) 



Remarkably, the charge and spin degrees of freedom are 
almost disentangled in the action S. The latter involves 
only the exchange interaction Hs- The remnant traces of 
the charging interaction He are encoded in the variable 
4>o, leading to a small imaginary shift of the chemical 
potential. 



III. WEI-NORMAN-KOLOKOLOV 
TRANSFORMATION 

The evaluation of the Green's function Q a {Ti2) is more 
convenient to perform in the Hamiltonian formalism. 
Then it can be written as 

g ,s_^_\-fcot<Tia 2 {-iT,-iT + iP), T > 0, 

Z I K. aaicr2 (-it - 0, -it), t < 0. 



Here Z = Trexp(— and we introduce 



K-acri cr 2 (^+ : ) Tl e 



(19) 



(20) 



Next, using the following set of transformations for the 
evolution operator (we recall that Hq and Hs commute) , 
we write 



N 



a it6 n s a /N 



hm T\e^ tJS ' N 

JV->oo J -- L 
n=l 

N 

n=l J a 

= [ V[G]e ± ^So dt ' e2 \[Te l K dt ' es % (21) 

^ a 

and obtain 



p=± 



ri4 +) <c 



(22) 



Here we have introduced the bosonic fields p , p = ±, 
and 



= Tcxp ( i / cft'0 p s Q 



(23) 



In Eq. ([2~Tj) we have missed the correct normaliza- 
tion factor (depending on J) because in the Hubbard- 
Stratonovich decoupling we were not pedantic enough 
concerning the normalization factor. For the computa- 
tion of the Green's function it is irrelevant due to the 
cancelation of normalization factors. For the partition 
function we restore it later on by comparison with known 
limiting cases. Note that while Ji is time independent, 
the factors Aa involve time ordering (T). This is due 
to the non-commutativity of the spin-operators s a . 

In order to overcome the intricacy of time-ordering we 
apply the WNK transformatio n 16 ^ 7 of variables in the 
functional integral (|22|) (sec Appendix [S] for details): 



Pp ^K P K p P . 



dl - ipOl 



— ipk p + p p n p — (n p ) 2 K p p 



(24) 



Here new variables p+ , k~^_ , k + correspond to 9 + whereas 
new variables p_ , kT , k~Z are introduced instead of 0_ . 

The WNK transformation recasts the time-ordered ex- 
ponent as a product of simple Abelian ones: 



exp 



exp 



is* / dt' Pp (t') 



x exp 



(25) 



Here = ± is^., and we employ the initial condition 
k p (0) = (the origin of this initial condition is discussed 
in Appendix • We stress that Eqs ((Ml) and (J2SJ are 
valid for a general spin operator. Originally, the field 
variables 6 p were real, but before the change of vari- 
ables (fM]) we have rotated the contour of integration in 
the complex plane. This procedure does not interfere 
with convergence of the Gaussian integrals. In order to 
preserve the number of field variables (three) we impose 
the following constraints on the otherwise arbitrary new 
complex variables: p p = —p* and k+ = (n p )* . The Jaco- 
bian of the Wei-Norman-Kolokolov transformation (|2~4"]) 
is given as (see Appendix El) 



j= n° x p(? f p dtp P (t) 

P =± V 1 J « 



(26) 



In terms of new variables the quantity /C ao - lCT2 (t + , t_) can 
be then rewritten as 



JC aoioa (* + ,t_) = J] / V[ Pp ,4 P ]e^fo Pdt ^ p ^ 
P =± J 

xe f/o p %(*) Ciff3 (t +i L) J] 7 (t+,i_), (27) 
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with C Q(Tl(T2 and B 1 given in terms of single-particle 
traces: 



B-y = tr 



-iM+^(+)( t+ ) e ^-y*-^(-)(t_) 



(28) 



Here h a = e aM n atCr . The expression for Z can be ob- 
tained from Eq. (|27[) by the substitution of B a for C aaia2 : 

, e f / tp dt Pj) (t) Y[B 7 (t + ,t-). (29) 



X( 



Simplifying expression (|25[) with the help of identity 
(s?J) 2 = (valid for spin 1/2) and evaluating the single- 
particle traces, we find the following result: 



I _|_ e -2ie 7 (t+-t_) 



2e -<Mt + -t-) C0S [l£ | 





p= 






exp 


^2 y 








/ d*«?.p 








dtp p (t) 



ip f dt' P-p(t') 

= o 



.(30) 



Here the presence of Zeeman splitting is taken into ac- 
count by means of the variables (b = g^sB) 



P P (t) = P P (t) - pb, k p ±p (t) = n p ±p (t) 



(31) 



The evaluation of the single-particle traces yields the 
following non-trivial matrix structure of C a in the spin 
space: 



p=± 

.4. = e -+ [ie <£ "<**+ei /o + dt ~P+M j + dt'k+(t') 



p c — 2ie Q t 



(i 



x e i/o' drp_(r) + e ie t+ £+(f') e f /o + dtp+(t) 



P =± 



x e 



lp / tp <itp p (t) f" difa-P^e-iPlo drp p (.r)] ( 32 ) 



We emphasize that the WNK transformation explicitly 
breaks the symmetry S z — > — 5 Z together with b — > — b. 



For example, C a \^. cannot be obtained from C a ^\ by re- 
versing sign of the magnetic field b. We shall see below 
how the symmetry restores. 



IV. EXACT EXPRESSION FOR THE 
PARTITION FUNCTION 

The partition function Z is given by Eq. (|2"9"]l . We 
start from integration over the fields n v p and K^L p . The 
expression (f3T))) for 2? 7 is bilinear form of the fields nP p and 
n p _ p . By using the following identity 



dz e 



2iti z 2 



(33) 



we rewrite Eq. (|2"9")l as 

Z = II/ P[p P ,^ p ]e-^/o p *(^-^ K --2^(t)) 



p=±' 



n 



"7 

2ttz2 

7 Kl=l 



exp ( — z*yB~/(t+, t- 



(34) 



Then the functional integral over the fields ftp p becomes 
Gaussian. As shown in Appendix [Bj due to the specific 
form of the initial conditions, they have simple dynamics 
and can be integrated out exactly. The result is 



z=J 2 (n fv[ Pp ] 
\p=± j 

n / ^ 

7 Kl=i 



e 4 



^J7 dtp 2 p + !?C dtp p (t) 



-w — 2v cos 



i E /dtp p (t) 

p=± 



exp 



-y-iJy II e 2P/o 

n/ 



!f S„" dtp p 



X U 



,V P eT('H-) f" dte -ip^dt' Pp (t') 
p=± •'o 



(35) 



where 



]Tz 7 Jl + e -2«7(*+-*-) 



(36) 



To transform Eq. (|35l) to more standard form, let us 
introduce new variables 



c P (t) = i P / dt'ft,^) + e P (o). 

Jo 



(37) 
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Here £ p (0) is an arbitrary constant. Then the partition 
function Z can be written as 




n / v m 

p=± j 



x cxp ■ 



-2ucosh 9^p(*p) - £p 



Vp^ 2 ' 



(0) - ibt p ] 



(38) 

The functional integral over fields £ p in Eq. ([58)1 is of the 
Feynman-Kac type. The quantity 



\<?=± 



3 ?P 



(39) 



plays a role of Lagrangian. It is convenient to perform a 
shift of variables £ p and introduce new variables 

i P (t) = - \ J2 &(*«) + M€«(0) + ipbqtq] ■ (40) 

Then the expression for the partition function Z acquires 
exactly the same form as given in Eq. (|38[) with the fol- 
lowing substitutions: £ p (t p ) — > £ p (t p ), £ p (0) — > £p(0), 
C p — >• £ p , where 



£ - — £ 2 - 



(41) 



We mention that the new variables £ p are, in fact, inde- 
pendent of the values of £ p (0). They obey the following 
constraints 



'%2ij > (t p )+2hi4yv = 0, 



P =± 



(42) 
(43) 



In what follows we shall omit the tilde signs on the vari- 
ables £ p . At this point it is convenient to express Eq. ([55)1 
in terms of matrix elements for the one-dimensional 
quantum mechanics with the Hamiltonian 

f) 2 J 

H.j = -J^ + -e-t. (44) 



Then the partition function Z becomes 

" J 2 dy 



n 

V 7 i* T l 




- w ( n / 

x e -2„cosh[(f + -f_)/2] ^ + | e -iff ji+ |£\ | e iffjt- 



(45) 



The Hamiltonian Hj of the one-dimensional quantum 
mechanics is exactly solvable. Its eigenfunctions are 
spanned by modified Bcssel functions K 2 i U '- 



HO = - v / ^sinh(2^)X 2w (e-«/ 2 ), 



(46) 



where v is a real parameter. The corresponding eigenval- 
ues of Hj are equal to Jv 2 : Hj\v) = Jv 2 \v). 

Next we perform integration over y in Eq. (|45p . Then 
with the help of the follo wing identity (see formula 
6.794.11 on page 794 of Ref. [If) 



dv v sinh(27w) K 2iv (e~^ /2 )K 2w (e~*- /2 )K 2w (2v) 



exp 



16 V 4v 
we integrate over £ p and obtain 



2t> cosh ■ 



(47) 



/idz I e~ w f°° 
\ - — / dv v sinh(27r^) 
2nz 2 v J 

ffao,(2t;) ^11 / <e-^/ 2 ^(e-4/2)j 




Vp=± / 



-iji/ 2 (t + -t_) 



(48) 



Next, u sing the identity (see formula 6.521.3 on page 658 
of Ref. Haf) 



dxxK u {ax)K u {bx)= ^ ah \ ( ° \} , (49) 
2sm(7w)(a^ — b 2 ) 



we perform integration over £' + and . With the help 
of the well-known integral representation of the modified 
Bessel function 



K 2w {2d) 



1 



dh e 



-2d cosh h+2iuh 



(50) 



we integrate over the variable v. Finally, integration over 
Zj can be performed, and we find 



Z = 



Hp 



-pb 2 /AJ 



x e 



-h 2 /pj 



dh sinhm) , , 'J 
v y sinh(/36/2) 

-Pe~,—ha\ 



(51) 



During the set of transformations we omitted normal- 
ization factors which depend on the parameter J. In 
order to restore them, one can compute the partition 
function Z for a single and two-level cases. Then one 
finds that the following transformation is necessary 



J 2 ' 



(52) 
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Hence, we obtain the following result for the partition 
function corresponding to the Hamiltonian %: 



1 



-/3(6 2 + J 2 )/4J 



x e 



-h 2 /pj 



IK 1 



' — oo 



ah sinh ft — -— — - 
v ; sinh(/%/2) 

(53) 



7,(7 



With the help of Eq. (fl3jl , the grand canonical parti- 
tion function for the full Hamiltonian ^ can be written 
as 



Z = 



VP 



-/3(b 2 + J 2 )/4J 



E< 



-/3E c (n-N y 



7T//3 



x e 



»»» r & sinh(/ l )^^e-'' 2 /^^ 
i-oo 7 sinh(/?6/2) 



-^n o ( M -i0 o +/i<T//3) 



2- 



(54) 



To integrate over the variables <po and h we can use the 
following identity for the grand partition function of free 
electrons 



no 



E z 

N=0 



N e 



(55) 



where the canonical partition function of N noninteract- 
ing spinless electrons is given by Darwin-Fowler integral: 



d6 . 

ZN =lo ^ 



-i6N 



](l + e ie -^) 



(56) 



Hence we find another representation of the grand canon- 
ical partition function for the Hamiltonian 



z= E 



bllm 2 7 7 -pEdn-Nof+pfin 



sinh ■ 



x e 



(57) 



Here n^(n^) represents the number of spin-up (spin- 
down) electrons, the total number of electrons n = 
nf + ni, and m — (n-f — nj,)/2. Note that for m ^ 
(m < 0) the total spin S = m (S = —m — 1), respec- 
tively. Different terms in Eq. (|57[) have clear physical 
meaning. The quantity E c {n — Nq) 2 — Jm{m + 1) is the 
interaction energy of the state with and electrons. 
The factors Z n ^ and Z ni take into account the contribu- 
tions from the single-particle energies. The 6-dependent 
factor 



sinh 



(96(27)1+1) 



sinh 



fib 
2 



5] cxp(/?6S 2 ) (58) 



S z — — m 



represents the partition function for spin S = m in the 
presence of Zeeman splitting. Finally, we mention that 
Eq. ((57)) coincides with the result obtained in Refs. 0(2^| 
by other approach. 



V. EXACT EXPRESSION FOR THE 
TUNNELING DENSITY OF STATES 

To derive an expression for the TDOS we begin from 
evaluation of the correlation function /C acrio - 2 which is 
given by Eq. (f2"T)l . As in the previous section, we start 
from integration over the fields K p p . The quantities C a -\i 
and Cq^t are nrs * order in n^ p . Therefore, the corre- 
lation functions K, a ^i and K, a ]/[ vanish after integration 
over the fields K p p . There is difference between C a -j~j- and 
Can ■ The former is independent of the fields k^ p whereas 
the latter does. Such an asymmetry is due to our choice 
in parameterization of the time ordered exponents (see 
Eq. (|25p ). In what follows, we shall evaluate the correla- 
tion function K, a n- Then integration over the fields k^ p 
can be done in exactly the same way as in the previous 
section for Z, since the quantity C a -\-\ is independent of 
the variables k^ p . We thus obtain 



n 

V* a |*yl = 



2^2 



-w a —2v a cos 



I E fdtp p (t) 

p = ± 



: / dye v exp 
'o 



x e 



2ie a t + gie„tp e -f |J P dtp v (t) ^ 



(59) 



where 



P =± 



> Ja = ^z 7 (l + e- 2 ^ ( *+-*- ) ) 



7^q 



u « = E 



y(t+-t-) 



(60) 



As in the previous section we perform a transformation 
of variables from p p to £ p and write the result in the 
Hamiltonian formalism (omitting tilde signs): 

"°° Pdy 



n f 



2-kz 2 ) J iyv a 



-y—w ct —2ie a t-i r . 



X <5 ^E ^ + 21n4y WQ ^J <5 ^ p% + ibt p ]j 

x e -2 Ua cosh[((; + -f_)/2] ^ + | e -ii?J«+ | e iffji_ |^_^ 

(61) 
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Next we perform integration over y in Eq. (|6"T|) . With the 
help of the identity (|47[) we obtain the following result 



K.n— in / ^ £■ 







#2i*O a ) ^ 

p=± 



x (He« /2 h)Q 



(62) 



where 



Qi/i^i (-z) = — o [ w i sinh(27rv) sinh(27ri/i)] 



1/2 



x / d V r ] 'K 2lv (T 1 /z)K 2wi (zT 1 ). (63) 
Jo 

Usin g th e identity (see formula 6.576.4 on page 676 of 
Ref. 11), we find 



ct,(t' — ± 



r ( - + iau + ia'v\ 



'3 3 
x 2-F1 ( - + iv + ivi, - + iv — iv\i 3; 1 — z~ 



(64) 



where T(x) and 2 F\{a 1 b, c; x) stand for the Gamma and 
hypcrgcomctric functions, respectively. With the help of 
the following relation between the modified Bessel func- 
tions 



we evaluate the matrix element as 



{v~i/2\ Vl ) | {v + i/2\ Vl ) 



(65) 



(66) 



Combining Eqs we obtain 



-2ie Qt t+ /3fa/2. 



d/i sinh(/i) 



2V7T? 

x e tpJt " /4 e zpi ^ t "W(2h + i P Jt p , P l3b/2,l3j) 
P =± 

x II II i 1 + e"^""' 1 ) • (67) 

Here we use the fact that t + — i„ = — i/3. The function 
W is defined as 



4 sinh y sinh y 



oJlrz ( x — 2ay 
erh ' 



2^~z 



4e y exp 



£ - 2y 
2^1 



(68) 



where erfi(z) = (2/y/n) f* dt exp(t 2 ) is the error function 
of an imaginary argument. The expression for IC a n can 
be obtained from Eq. (|67p by transforming b — > —b. 
Finally, from Eqs JT^]) and (Q3]) for t > we obtain 



G «tt( T ) 



1 \/? r _j. 



/3B c (n-JV )2 



x e 



-B c (2n-2AT + 1)7 



^0 e i/30 o n e -( £clT - AJ )7 



7T//3 



<i/isinh(/i) J^J - 



2tt 

/3U (tJ.-i<t>o+ha / 0) 



a=± 



_|_ e -P(e a -fj,+i<p )+ha 



x |e J(2r -«/ 4 W(2/i + Jr,j3b/2,J3J) 

+ e -P(^-»+i<t>«) w (2h + J{P - t), -fib/2, PJ) | . 

(69) 

The integration over (j)o in Eq. (|69p can be performed 
with the help of the identity 

oo 

I (l + e-> 3 ^- p A = Z N {e a )e^ N , (70) 



N=0 



where Ziq{e a ) is the canonical partition function of a sys- 
tem of N noninteracting spinless electrons under the con- 
straint that level a is not occupied. It is given by the 
Darwin-Fowler integral 



2 71 



de 

2V 



— iQN 



(71) 



7^q 



and obeys the following identity: 

Z N = Z N {e a ) + e-^Z N ^{e a ). (72) 



Performing integration over 0o and h in Eq. (|69[) . we 
obtain (for r > 0) 



«ttW = -2z E 



-0_E c (ra-Aro) 2 +/3/in+^Jm(m+l)] 



X e -[ c «t-M+Bo(2ri-2W + l) + J(m+l/4)]r| e /3()/2 

x T(/3b, 2m + 1) Z nf (e a )Z rH - Z„ t+1 Z ni _i(e a ) 
- T(-(3b, -2m) Z„ T Z IU (e Q ) - Z„ T (e Q )Z, H j, (73) 
where 

T( , xl . e( ^ lW2 sinh(xz/2) 
V ' ' sinh(z/2) a;sinh 2 (z/2)' V ' 

We notice that Y(z — > 0, x) = x — 1. The expression for 
G a jj (t) can be found from Eq. (|73| by reversing the sign 
of the magnetic field b. 
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Employing the general expression for the TDOS^l 



■ — cosh — 

7T 2 



die 1 



(75) 

finally, we find the following exact expression for the tun- 
neling density of states for the UH @: 



/3fe(2m+l) 



v a {e) = 



7 7 p 



0Jm(m+l) 




1 - S_ m _! (a(m + l)pb) 
S(s - e aa + /i - E c (2n - 2N + 1) + J(m + 3/4)) 

a ) ^n-t- 



(2m + 1)Z„ T 



1 + B m (amPb) 



S(e - e aa + fi- E c (2n - 2N + 1) - J(m + 1/4)) 



Here 

2m + 1 / 2m 
B m (x) = — coth 



2m 



2m 



(76) 



— coth-^- (77) 
2m 2m y ' 



denotes the Brillouin function. Equation ([76)) constitutes 
the main result of the present paper. It allows to com- 
pute the TDOS for a given realization of single-particle 
levels. As expected, according to Eq. (|To) . the TDOS 
represents a sum of delta-functions corresponding to all 
possible processes of tunneling of an electron with energy 
e and spin a into (or from) a single-particle level with en- 
ergy e aa . The factors Z n (e a )/Z n describe the probability 
that the single-particle level a is empty. 

By using the identity J2 a l z n - Z n (e a )]/Z n = n, one 
can check that the result (fTBf satisfies the sum rule: 



de 



1 + exp(e/T) 



= -oT- 



dlnZ TdlnZ 



db 



2 dfi 



(78) 



In the case b = J = and for spinless electrons the 
result (|76p coincides with the expression for the TDOS 
found in Ref. 



VI. STATIC SPIN SUSCEPTIBILITY: THE 
EFFECT OF LEVEL FLUCTUATIONS AND 
ZEEMAN SPLITTING 

In this section we consider the thermodynamics of the 
quantum dot at relatively low temperatures, S « T « 
E Cl i?Th- The quantity of main interest is the static spin 
susceptibility averaged over realizations of single-particle 



levels in the presence of Zeeman splitting. As is well es- 
tablished, its divergence indicates the Stoner instability. 
In general, the static spin susceptibility is defined as 



x(T,b)=T 



d 2 hxZ 
db 2 ' 



(79) 



In order to compute it one needs to perform integration 
over </>o and h in Eq. (f54"| . At temperatures T ^> 5, 
the integration over </>o can be performed in the saddle- 
point approximation! 1 ^ 12 Then, the partition function 
becomes 



with 



Zs 




Z — ZcZs-, 



-/3B c (n-AT ) 2 + / 3 Mo n-2 / 3n (A) 



(80) 



(81) 



dh sinh(/i)e 



-h 2 //3J 



sinh(/36/2) 



(82) 



Here jj, = fi + with fig being the solution of the saddle- 
point equation: 



N 



(83) 



This means that at T S the charging and the spin 
part of the problem are decoupled from each other. The 
partition function in the absence of the exchange interac- 
tion and magnetic field is given by Zq- The factor Zs in 
Eq. (|5D|) describes the effect of the exchange interaction 
and magnetic field: Zs — 1 at J = b = 0. Next, 



A = - 



dfi. 



r,2 



(84) 



stands for the inverse thermodynamic density of states 
at the Fermi level for a given realization of the single- 
particle spectrum. Since Zq is independent of the mag- 
netic field and exchange interaction it does not affect the 
spin susceptibility. Therefore, in what follows we will 
discuss Zs only. In the absence of exchange interaction 
InZs equals to (3 ^ (T [fio(A) — ^o(A + ba/2)] as it should. 

The function /3^ CT [r2o(A) ~ ^o(A + her /ft)] that ap- 
pears in Eq. (|82p is a random function of the variable h 
due to fluctuations in the single-particle density of states 
vq{E) = J2 a + p, — e a ). Provided h 2 < exp(j3jj,), we 
find 



/3^[f2 (A)-fioO« + W/3) 



(35 



V(h), (85) 



where 
V(h) = 



dE5vo(E) In 



sinh 2 (/i/2) 



cosh 2 (.E/2T) 



(86) 
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is a random function. Here 8v$(E) stands for the devi- 
ation of the single-particle density of states Vq (E) from 
its average value: 



1/5=1/ A. (87) 
With the help of Eq. (JSS) we rewrite Eq. (JS5]) as 



1 E(6/J,/3J*) 



exp 



V^TP sinh(/36/2) 

/OO 
d/i sinh(x/i)e 
-oo 



-/3 



J 2 + 6 2 



4./ 

h-h 2 / V -V(h) 



Here 



J. 



J<5 



(88) 



(89) 



denotes the renormalized exchange energy. We empha- 
size that near the Stoner instability, 6 — J -^i 6, the value 
of spin in the ground state is of the order of 6/ [2(6 — J)], 2 
and the renormalized exchange energy is much larger 
than the bare exchange, J* ^> J. In what follows, we 
consider this most interesting regime. 

Let us neglect V(h) in Eq. (|88[) for a moment. 
Then, the typical value of h in the integral of 
Eq. (|88[) is of the order of max{ v /y, y, yx}. There- 
fore, Eq. (|85j) is valid if the following condition holds: 
max{l, 1/^fjJJl, b/ J} <C exp(/3/i). It is satisfied for 
large enough values of the chemical potential. 

Although the single-particle density of states v (E) has 
non-Gaussian statistics, for max{|/i|, T} 3> S the function 
V(h) is a Gaussian random variable^ Its statistics is 
fully determined by the pair correlation function 



C(h 1 ,h 2 ) = V(h 1 )V(h 2 ). (90) 

The following exact relation holds (see Appendix |C|) : 

C{h u h 2 ) = L{h x + h 2 ) + L{h x - h 2 ) - 2L(h{) - 2L(h 2 ). 

(91) 

At T ^> S the behavior of the correlation function L 
depends strongly on the value of \h\: 



LW n*(3 \H\h\/2) 



\h\ « 1, 

c 2 - 3/2, 



\h\ > 1. 



(92) 



Here we introduce the parameter (3 such that (3 = 2 for 
the unitary Wigner-Dyson ensemble (class A) and (3 = 1 
for the orthogonal Wigner-Dyson ensemble (class AI). 
The numerical constants arc 



ci 



du I 1 u) coth uj — 1 
— \ 3 



1 duj 



sinh 2 u) 



f du> r , . f 

c 2 = — — t 1 - uj coth lj\ + 

Jo w Jl 




In spite of the fact that V(h) is a Gaussian random 
variable exact evaluation of In S(a;, y) for arbitrary values 



'. ~< 

! » 
i — ^ 












II a \Hb 


III 




\ 
\ 
\ 






X 
\ 
\ 
\ 
\ 








/ 


y=r/ x 



b/J 



FIG. 1: Different regions for behavior of the spin susceptibility 
in the plane of dimensional parameters b/J and J*/T. Note 
that in our analysis T 3> S. 



of x and y is a complicated problem (see e.g., Ref. [29||). 
In what follows, we compute lnS(a;, y) to the first order 
in the correlation function C. We expand expression ([55)) 
for S(x, y) to the second order in V and perform averging 
of lnS(x, y) with the help of Eq. ([9"Tj) . Then, we find 

In S(a;, y) = In S (.t, y) + i V / -^e"" 2 

x | cr coth + 1 i(y(l + xcj) + 2uy/y) 

- L(y(l + xa) + uv^y)] - L(u^) 
L(y{l + xa) + u^/2~y) 



2sinh 2 (a:y/2) 
- L(y + u^/2y) - L(yxa 

+ L(u^2y~) 



(95) 



where 



e v(i+. 2 )/4 sinh ^. 

2 



(96) 



Exact integration over u in Eq. (|95|) is complicated since 
only asymptotic expressions (|92|) for L arc known. To 
make further analytical progress it is useful to consider 
separately regions in which arguments of the functions L 
involved in Eq. are either large or small. This way 
one finds severa l regions s hown in Fig. [TJ in each of those 
the behavior of lnS(a;, y) is different. 
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A. Region I: J* max{l, b/J} < T 

In the region I, max{y , xy} -C 1, the arguments in 
all functions L in the right hand side of Eq. (|!l5j) are 
much smaller than unity. Therefore, in order to evaluate 
integrals over u we can use the asymptote of L(h) for 
\h\ < 1 (see Eq. fl§2J)). Then, we obtain 



lnS(x, y) = In S (a;, y) + 



Cl 



96/3tt 2 

+ 12xV(l + 2//2)-xV/6 



30y 2 + 12y 3 



(97) 



Hence, using Eqs (fT9")l and ([55]) . we find the average spin 
susceptibility in the region I (J* max{l, 6/ J} <C T): 



x(T ' 6) = 4{ 1+ 6^ + |^ 

Jfo 2 r 10c Jl_ 
120T 3 J 2 i + /3 T 2 



1 + A 
2T 



(98) 



Expanding Eq. (|100[) to the forth order in b and per- 
forming averaging over level fluctuations by means of the 
expression ~J™ = J™[1 + n(n + l)cJ 2 /(2/3T 2 )] we ob- 
tain Eq. (|98p . The expression (j9"8")l underlining the diver- 
gence at the Stoner instability point represents exchange- 
enhanced Pauli spin-susceptibility (J+/2J8) with small 
corrections depending on temperature and magnetic 
field. The corrections due to level fluctuations are small. 

It is worthwhile to mention that the spin susceptibility 
at b = has been studied by Kurland ct al.— In our 
notations, their result at T J* can be written as (see 
Eqs (4.8), (4.13b), (4.15) of Ref. 0): 



Xkaa(T,0) = A 



~6~Vf + CKAA T 



, (103) 



4 - 7T 



ckaa 



.(■ 



!ln2 



- 0.3712 



6 6/3 V 7T 2 
The result (|103[) of Ref. 0] disagreees with our result 



Here the numerical constant c = ci/2tt 2 rj 0.02. 

It is instructive to derive Eq. (|98p in a more transparent 
way. In region I the typical value of h in the integral in 
Eq. (|88[) is much smaller than unity. Therefore, we can 
expand the random function V(h) as a series in h. We 
thus find: 



V(k) 



< i. 



(99) 



Next, performing integration over h in Eq. (|88[) . we ob- 
tain 



B. Region II: S < T < J* and tcJ 

In the region II, y 3> 1 and a; 1, the arguments of 
all functions L in Eq. (|95p are typically much larger than 
unity. However, the behavior of some contributions to 
the integral over u in Eq. ((95]) are different for x 2 y 3> 1 
and x 2 y <C 1. Therefore, it is convenient to split region 
II into two regimes. 



J smh[Jb/(2JT)] j J _ j )(J 2 +b 2 )/iJ 2 T 



VI sinh(6/2T) 



(100) 



1. Region II a : S < T < J* and 6 2 /J 2 < T/J* 



where \jj = 1/A stands for the renormalized 

exchange interaction for a given realization of single- 
particle spectrum. All information about level statistics 
is contained in fluctuations of A. At T 3> 6 the fluctua- 
tions of A are small and Gaussian with (see Appendix [Cjl 



(A 



c S 4 

' (3T 2 



(101) 



To avoid the Stoner instability the renormalized ex- 
change energy J should be positive for a given realiza- 
tion of single-particle spectrum. Since we are interested 
in the regime J* 3> 5 > J, in order to fulfill the condition 
J > the level fluctuations should be such that 



(102) 



Using Eq. (jlOip . we find that Eq. (|102[) is equivalent to 
the condition T 3> Provided the latter is satisfied (i.e 
in the region I) the level fluctuations are small enough 
and cannot drive the system to be Stoner unstable. 



Let us first consider the regime of weak magnetic fields: 
x 2 <C 1/y <C 1 (region II a ). In this case, we perform 
integration over u in Eq. (|95p either by expansion of L 
to second order in u or with the help of the following 
asymptotic result at z 1: 



°° du _ u 2 

— e L(zu) 



2/3tt 2 



z 1 + 7 

ln 4 +C2 -^- 



(104) 



Here and below, 7 rj 0.577 . . . stands for the Euler's 
constant, and constant C2 is defined in Eq. (|94|) . We 
obtain 



lnS(x,y) = lnS (a;,y) 



1 



2/37T 2 



xy coth 



•fy 



9 2 

x y 



4sinh^(xy/2) 



(ln2y + 7 ) 



(105) 



Using Eqs (|88|) and (|105|) . we find the average spin 
susceptibility (defined in Eq. (|T9"|) ) in the region II a 
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(b 2 /J 2 < T/J+ « 1): 



2. Region II b : i « T « J 4 and J+/T < fe 2 /-^ < 1 



1 - 



Jib 2 



AJ 2 T 2 sinri (J*6/2JT) 



t a 2 



J*b J*b 

coth 

JT 2TJ 



2(3ir 2 db 2 

J 2 & 2 (ln(2J*/T) + 7 ) 
AJ 2 T 2 smb 2 (J ir b/2JT) 



(106) 



In the limit of very weak magnetic fields b <C JT/J* 
we can neglect tiny dependence of the average spin sus- 
ceptibility on magnetic field. Then, from Eq. (|106[) wc 
obtain the following result for the zero-field average spin 
susceptibility: 



X(T,0) 



12T 



1 

W 2 



In — + 7 + 2 



(107) 



The result (|107|) is valid provided the expansion of 
In y) in powers of the correlation function (|90|) is jus- 
tified. As one can demonstrate (see Appendix|D|), the lat- 
ter is controlled by the small parameter J*/{(3-k 2 T) <C 1. 
Therefore, strictly speaking, the result (|107[) holds at 
1 <C J+/T -c (3tt 2 . However, the more detailed anal- 
ysis of x(T, 0) presented in Appendix [Dj together with 
the qualitative arguments of Sec. IVIII below, allow us to 
speculate that the result (|107p has a much broader range 
of applicability. 

We expect that the zero-field average spin susceptibil- 
ity at J* T 3> S can be written as 



X(T,0) 



{J* I J? 
12T 



1 

W 2 



hi 



J * 



•J jf 



V/37T 



/3tt 2 T 



(108) 

where f(x) is constant of the order unity in both limiting 
cases of small (x <C 1) and large (x ^s> 1) values of its 
argument. In particular, Eq. (|107[) implies that f(x) — 
7 + 2 + In 2 H at x < 1. 

Therefore, near the Stoner instability, 5 — J <C 5, 
there is ehnancement of the average spin susceptibility 
(Eq. (|107[) ) due to fluctuations of single-particle levels in 
a wide temperature range 5 <C T <C J*. 

In the regime of larger magnetic fields JT/J* <C b <C 
sfjT/Ji from Eq. (|106p we find the following result for 
the average spin susceptibility: 



x(T,b) = ^ 



Jib 2 

J2 T 2 



1 J 2 J 2 
2/3n 2 J 2 T 2 



-J*b/JT 



(109) 



The spin susceptibility is suppressed as compared to the 
zero- field result (|107[) . It is given mainly by T/b 2 term 
whereas the level fluctuations contribute to the terms 
which are exponentially small, oc exp(— J+bj JT). 



Let us now consider the regime of intermediate mag- 
netic fields: 1/y <C x 2 <C 1 (region lib). Performing 



integration over u in Eq. (|95|) , as in the previous section, 
either by expansion of L to the second order in u or using 
the asymptotic result (| 104f) we obtain 



\nE(x,y) = \nS (x,y) + 



2 2 

x y 

W 2 



111 



x--\e 



(110) 



With the help of Eqs flSS) and ([LTU]) . we find that the 
average spin susceptibility in the region II;, (J/J* <C 
6 2 /J 2 <C 1) is given as 



2J* 
(3ir 2 T 



b 
J 



(111) 



We note that a small magnetic field b ~ (J/J*)T < J,T 
destroys the phenomenon of mesoscopic Stoner instabil- 
ity, i.e. the existence of the ground state of a quantum 
dot with a finite non-zero spin. It renders the spin- 
susceptibility of the Fermi-liquid type. Since J*/ J is a 
factor responsible for the enhancement of a g-factor in 
the Fermi-liquid, the energy scale J+bj J determines the 
effective Zeeman splitting. As usual, comparison of the 
latter with temperature allows us to distinguish between 
weak and strong magnetic fields. 



C. Region III: i>J*/J and J < b 

Similarly to the region II, in the region III (y ^> 1/x 
and x ^> 1) the arguments of all functions L in Eq. (|9l))) 
are typically much larger than unity. However, in the 
region III for evaluation of the integral in Eq. (|9"5|) one 
can expand all functions L to the second order in u ex- 
cept L(uy/2y). The latter can be integrated by means of 
Eq. ([TD3J) . Then, we obtain 

In E(x, y) = hi Ho (a, y) + (i n ^ + c 2 ) . (112) 

Hence, using Eqs ([581 . we find the average spin suscepti- 
bility (defined in Eq. ([79])) in the region III (TV J* < 6/ J 
and J < 6): 



(113) 



We recall that in the course of the derivation of Eq. (jl 13|) 
we have neglected terms exponentially small in b. The 
spin susceptibility (|113p is exchange-enhanced Pauli sus- 
ceptibility, with correction due to fluctuations of single- 
particle levels. The latter is small in the regime consid- 
ered. It suggests that taking the lowest order expansion 
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in V(h) for evaluating lnS(a;, y) is well-justified in the 
whole region III. We mention that our result (|113[) for the 
spin susceptibility is an extention of the zero-temperature 
result of Ref. [1] with account for the effect of level fluc- 
tuations. 

It is instructive to summarize here the results of the 
above analysis of the average spin susceptibility. In re- 
gion I (see Fig. [T]) , the mesoscopic Stoner instability is 
manifested through small temperature and magnetic field 
dependent corrections to the Fermi-liquid result. The ef- 
fect of level statistics is weak. In regions lib and III, 
the mesoscopic Stoner instability is suppressed by mag- 
netic field. The corrections to the Fermi-liquid result 
due to fluctuations of the single-particle levels are small. 
In region II a , the average spin susceptibility behaves in 
accordance with the Curie law. The latter is a mani- 
festation of the mesoscopic Stoner instability at non-zero 
temperature. Fluctuations of single-particle levels lead to 
logarithmic-in-tcmpcrature corrections to the Curie term 
in the average spin susceptibility. 

Here we do not study fluctuations of spin susceptibil- 
ity due to statistical fluctuations of single-particle levels. 
However, we expect that in all regions except region II a 
they are small. In region II a , we estimate (x 2 — (x) 2 )/(x) 2 
to be of the order of (In J+/T)/ ((3tt 2 ). This result indi- 
cates that in region II Q fluctuations of the spin suscep- 
tibility can be large and, therefore, it is challenging to 
study the whole distribution function for the spin sus- 
ceptibility. 



VII. SEMI-QUALITATIVE ANALYSIS OF THE 
EFFECT OF THE LEVEL FLUCTUATIONS ON 
THE SPIN-SUSCEPTIBILITY AT T = 

A. Spin susceptibility at zero magnetic field 



The origin of the logarithm in the result (|107|) has sim- 
ple physical explanation. Although the result (|107|) was 
derived for T 5, let us consider the limiting case of 
vanishing temperature, T = 0. Then, the transition be- 
tween the ground states with the spin 5 and 5+1 occurs 
if 



E s+1 -E s = J(2S + 2), 



(114) 



where Es stands for the single-particle contribution to 
the ground state energy ('kinetic energy') for a given re- 
alization of single-particle levels. It can be estimated as 

-E s = 5(25+1) + AE 2S . (115) 

Here AE 2 s is a fluctuation of the energy strip in which 
there are 25 levels in average. We can estimate it as 



AE 2 s = SAn 2 s, 



(116) 
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FIG. 2: (Color online) Comparison between the average spin 
susceptibility at zero temperature and zero magnetic field es- 
timated from the theoretical result (|107[l with T = 5 (solid 
curve for (3 = 2 and dashed curve for (3=1) and results of 
numerical simulations of Ref. @] (circles for (3 = 2 and squares 
for (3 = 1). 



find 



2(5 - J) 



1 - An 2S 



(117) 



Hence, the average spin susceptibility at low tempera- 
tures can be esimated as 



X(T,0) 



5(5 + 1) 



J* 



l + (An, VJ )2 (118) 



3T 12TJ 2 
From random matrix theory it is well-known tha1j=£ 

2 



(An 2S ) 2 



/3tt 2 



In 25 + const 



(119) 



Hence we find the following estimate for the spin suscep- 
tibility at T < 5: 

J 2 r 2 / J* \- 
1 + ——-^ hi — + const 
L (3n z V J 1 - 



X(T,0) 



12TJ 2 



(120) 



The estimate f|l 20[) derived from qualitative arguments 
for T <C 5 resembles the result (|107[) which is valid at 
higher temperatures T 3> 5. There is a discrepancy 
in factor of 2 in front of the logarithm between esti- 
mate (|120p and the rigorous result f|10T[) . The reason for 
this discrepancy is the following. The qualitative argu- 
mens presented above correspond to evaluating the inte- 
gral over h in Eq. Q88JI in the saddle-point approximation. 
In such a procedure one misses the contributions of the 
type L(uy/2y) in Eq. (pj)) . 

In Fig. [5] we present the comparison between the av- 
erage spin susceptibility at zero temperature and zero 
magnetic field estimated from Eq. (|107p with T ~ J and 
results of numerical simulations of Ref. Q . 



B. Spin susceptibility at a strong magnetic field, 



where An 2 s is a fluctuation of the number of levels in 
the energy strip AE 2 s- Near the Stoner instability, we 



The result (|113[) can be illustrated by simple physical 
argumentation. As above, we consider the limiting case 
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of vanishing temperature, T = 0. Then, the difference 
between the single-particle contributions to the ground 
state energy in Eq. (|114|) can be estimated as 



E s+ i -E s = 8{2S + 1) - b + AE 2S . 
Near the Stoner instability, we find 



S = 



1 



2(8 - J) 



b - 8 An 2 s 



(121) 



(122) 



Hence, the average spin susceptibility at low tempera- 
tures can be estimated as 



x(T,b) 



dS_ 
db 



2J5 



1 



Jl d 2 (An z ) 2 



28 2 dz 2 

z = J t b/JS 

(123) 

Using Eq. (|119[) . one exactly recovers the result (|113D . 



Since our result (|113[) can be obtained from zero- 
temperature considerations presented above, we expect 
Eq. (|TTg|) to be valid at T = provided b > J, 8. 



VIII. TUNNELING DENSITY OF STATES IN 
MAGNETIC FIELD IN THE ABSENCE OF 
LEVEL FLUCTUATIONS 

In this section we analyze the TDOS at sufficiently low 
temperatures 8 «T <C E c , E^. As it was demonstrated 
in the previous section the effect of exchange interaction 
is most pronounced in the vicinity of the Stoner instabil- 
ity, 8 — J <C 8. Therefore, we shall consider the regime 
J* 3> J below. As was indicated in our above analysis 
of the spin susceptibility, here, too, there are the same 
three regions with different dependence of the TDOS on 
temperature and magnetic field (see Fig.[l|. The fluctu- 
ations of single-particle levels are important in region II a 
only. Most of our discussion below excludes the effect of 
level fluctuations. We eventually present semi-qualitative 
arguments to account for the effect of the latter. 

We start from rewritting Eq. as 



1 VP 



x e 



Jr/4 p (/3-r)&/2 p -l3E a (n-N 
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- E o (2n-2N + l)T+irj> 



?° Wo" 



7T//3 



2 71 



dh 



.<7=± 



- pn {p-i<j> +ho 1 13) 



G (t,/j,- ifio + h/P) 



x < e 



J(2T-f3)/A(l-T/p)h 



pW(2ph + Jt, pb/2, p J) 



P =± 



e -f3b/2 e -hT/0 J2 P W(2ph + J(P - r), -pb/2, PJ) 



P =± 



(124) 



Here the Green's function of non-interacting electrons in 
imaginary time is given by 



G (t,h) = -[ dE 



ME) 



2cosh(/3 J B/2) 



(125) 



Provided the condition 5 <C T <C \i holds, the Green's 
function (|125[) can be simplifcd to 

G (r, (i) = -(ttT/S)/ sin(^Tr). (126) 

Integration over 4>o and h in Eq. (|124p can be performed 
in the same way as it was done in the previous sec- 
tion for derivation of the spin susceptibility. We re- 
mind that integration over h was done under assumption 
that PJ+ max{l, 1/V/3J*, b/J} <C exp(/3/2). Then, using 
Eq. ([50]), wc obtain 



G n (r) 



tt/(P6) 



-0E c (n-N o +T//3) 2 



sin(7TT//3) e-^^C""^) 2 

nGZ 



,E c t(t/P-1) 



X<M + 



l^pje 



8yQ* cosh(/36/4) sinh(J*/36/2J) 
Y a \etft(^/p~Z(ab - J)/2J) 



-erft(yfpX(ab + J(1-2t/P))/2J) j. (127) 

The Green's function G^(t) determines the TDOS in 
accordance with Eq. (|75l) . Performing integration over t 
with the help of the following identity: 



dt 



cosh(7rf) cosh(z/2)' 



(128) 



we find 
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crfi 



2 J 



f F ( P s-2 P n-p + — 



-¥(p{pe - 2 P n- p + a-6/2), sb/2J, y/p£\ 



E- 



-PE c (n-N„) 2 



LriG2 



(129) 



Here vq = 1/8 denotes the average density of states of 
non-interacting electrons for single spin projection, f2^ = 
E c (n - N +p/2), f F (e) = [1 + cxp(e/T)]- 1 stands for 
the Fermi function, and the function 



F(x,y,z) = U-*/ 2 



dt 



cosh(7rt) 



erfi(z(y — it)) . 



(130) 
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FIG. 3: The tunneling density of states in the Coulomb valley. 
The solid (dashed) line corresponds to J/8 = 0.92, 8/T = 
0.35, and J*/T = 3.95 {J/5 = 0.92, S/T = 0.95, and J*/T = 
10.70). The inset depicts the nonmonotonic behavior. 
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FIG. 4: The tunneling density of states at the Coulomb peak. 
The parameters are the same as in Fig|3] The insets depict 
the nonmonotonic behavior. 



Equation (|129[) describes dependence of the tunneling 
density of states on energy, temperature and magnetic 
field. It is valid under the following conditions: 5 <C 
T <C /i and (3 J* max{l, l/y//3J+, b/ J} <C cxp(/3/2), and 
near the Stoncr instability, J* 3> J. We remind that 
fluctuations of the single-particle levels arc not taken into 
account in Eq. (|129|) . 

As follows from Eq. (1129[) . at small magnetic fields, 
b < J, the dependence of the TDOS on magnetic field is 
very weak. The TDOS for b < J almost coincides with 
the result for 6 = 0. Therefore, below we consider in 
detail the cases of zero (b = 0) and large magnetic fields 
(b > J) only. 



A. TDOS in zero magnetic field, 6 = 

Expanding the expression in square brackets in the 
right hand side of Eq. (|129[) to first order on b, we find the 
following result for the TDOS in zero magnetic field 



^(g) 



E - 



-/3E c (n-N ) 2 



l + —)f F (pe-2 P n-n 



J ( pe-2pW n 
2J/1 J* 



NoV 



(131) 



Here the function F(x, y) is defined as 



_F(x,y) = -e-^V*/ 2 



dt ■ 

x cosh(Trf) 



(132) 



Using the following asymptotic expression at y <C 1 for 
T(x,y): 



y 



2cosfT (yx/2) 



(133) 



we obtain the TDOS at T > J*: 



ty(g) 



E < 



-/3E c (n-N ) 



'f F {pe-2p£l-P){l 



4cosh"(|e-paT p ) 




(134) 



It is instructive to compare (|1 34[) with the result of 
Ref. [l3[ for the TDOS in the absence of exchange in- 
teraction. As expected, in this high temperature regime, 
the exchange interaction affects the tunneling density of 
states only slightly. It is worthwhile to mention that the 
correction to the tunneling density of states due to ex- 
change interaction is of the order of J/T rather than what 
one may expect from above results on the spin suscepti- 
bility (see Eq. ([55])). namely, J*/T. 

In the regime of intermediate temperatures 5cr« 
J*, it is convenient to use the following result for the 
behavior of T(x, y) in the limit y ^> 1 (see Appendix IB")) : 



F( x > y) = 2 sgn ( cos 



1 - erf 



7r.r 



_l_ e %( x ~\ x \) 



(-l) m e~ ylxlm+vm{m+1) 6(\x\ -2m- 1). 

(135) 



m>0 



Here 9(x) is the Heaviside step function (0(0) = 0), and 
the error function erf(z) = (2/y / 7r) J Q exp(—t 2 )dt. 

As x is varied for a fixed y, Eq. (|135[) suggests that 
J^(x, y) exhibits damped oscillations with a period 4 
(equivalent to an energy scale 4 J*). However, it is not 
the case. At y 3> 1 the function J-(x, y) is monotonous 
and close to the function 1/(1 + cxp(— y(x — 1))). The 
linear combination of two Fermi functions (standard one 
and one shifted in energy on J*) in Eq. (|131[) leads 
to the appearance of a maximum in the TDOS. The 
height of the maximum can be approximately estimated 
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FIG. 5: The total tunneling density of states v(e) — ^t( e ) + 
vi(e) in the Coulomb valley for the magnetic field b = 
(solid line), b — 2.75 J (dashed line) and b — 3. 25 J (dotted 
line). The other parameters are J/8 — 0.92, 5/T = 0.95, and 
J*/T = 10.7. The inset depicts the nonmonotonic behavior. 
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FIG. 6: (Color online) Tunneling of an electron with spin up 
(left) and spin down (right) into a quantum dot with finite 
value of spin in the ground state (see text). 



as [^ (J (e)/i'o]max — 1 ~ J/ J*. This additional structure 
in the TDOS reflects enhanced electron correlations due 
to the exchange interaction. 

In the case of temperatures T <§C E c we illustrate the 
non-monotonic behavior in the TDOS due to exchange 
interaction in Figs. [3] and |4] for the Coulomb valley (Nq is 
integer) and for the Coulomb peak (Nq is half- integer) . 



B. TDOS at a strong magnetic field, b ^> J 

Now let us consider the case of a strong magnetic field, 
t > J. As above we consider the regime S <C T <C 
J* near the Stoner instability, J <C J*, in which the 
most interesting behavior of the TDOS takes place. To 
simplify the expression (|129p for the tunneling density of 
states, we use the following asymptotic result 



¥(x,y,z) 



1 



^y-'/^Tpy-x/z^z 2 ), (136) 



which is valid for z 1 and y 1. It is worthwhile 
to mention that Eq. (|136l) works well already for y > 2. 
With the help of Eq. (JT35J!, we obtain from Eq. ([T2"!)| the 
following expression for the TDOS 



Me) 



-/3E c (n-N y 
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2 J* (6 — J) cosh(^6/4) 
- JT((p £ - 2pOP - ,hb/J)/.h, (3 J, 
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-/3B c (n-Ar ) J 



(137) 



As follows from Eq. (|137[) . the non-monotonic behavior 
of the TDOS due to exchange interaction survives in the 



presence of magnetic field. In Fig. \E\ the dependence of 
the total TDOS v(e) = v^(e)+v^(e) on energy at different 
magnetic fields is shown for the Coulomb valley. The role 
of magnetic field is of two kinds. At first, it suppresses the 
height of the maximum: [^(e)/2i/o] ma x — 1 ~ J5/(J*b). 
Secondly, the width of the maximum increases linearly 
(~ J*b/J) with the magnetic field. Finally, we note that 
the difference v^{e) — v\,{e) is small. 



IX. THE EFFECT OF LEVEL FLUCTUATIONS 
ON THE TDOS IN ZERO AND STRONG 
MAGNETIC FIELDS 

In order to develop qualitative explanation of the max- 
imum in the TDOS and in order to qualitatively under- 
stand the effect of level fluctuations, it is instructive to 
consider a Coulomb valley at T = 0. We start from the 
case of zero magnetic field. Let us consider tunneling of 
an electron with spin up to the ground state (with spin 
S) of a quantum dot in a Coulomb valley (see Fig. [B]) . 
The tunneling is possible only if the electron energy e 
exceeds E\ = E s+1 / 2 — E$ — J(S + 3/4). Here we recall 
that Es denotes the single-particle contribution to the 
energy of the ground state with spin S. An electron with 
spin down can tunnel provided its energy is larger than 
£ 2 = E s _ 1/2 -E s + J(S +1/4). Therefore, only tunnel- 
ing of spin-up electrons is allowed in the energy interval 
Si < e < At very large electron energies the tunnel- 
ing is insensitive to the spin of electron. It is natural to 



denote the energy £ 3 = E s+1 / 2 — Es + J(S +1/4) as the 
characteristic energy above which there is no difference 
in tunneling for electrons with spin up and down. With 
the help of Eq. ([TTgJ) . one finds £ 2 - £ x w J - ^E 2S -\ 
and £3 — ^2 ~ 2JS* + AE^s-i- To sketch the tunnel- 
ing density of states we employ the sum rule ([75)1 . In 
the Coulomb valley and at zero temperature it renders 
the integral / dev a (e) independent of J. The TDOS at 
T = is shown schematically in Fig. [7] The relative 
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height of the maximum in the TDOS can be estimated 
as (£2-£0/(2(£ 3 -£ 2 )) ~ 1 /( 45 ')- Near the Stoncr insta- 
bility, J* ^> J, the relative height becomes of the order 
of J/(2J*) if one neglects the effect of level fluctuations. 
This estimate is in accordance with the result (|13ip de- 
rived for i « T « J 4 . In this temperature regime, the 
feature in the TDOS in the energy interval £ y < e < £ % 
(see Fig. [7]) is smeared, as expected, since £2 — £1 <C T. 

Our qualitative arguments can be justified by di- 
rect evaluation of the TDOS at zero temperature in a 
Coulomb valley. For the ground state of a quantum dot 
with spin 5 we find from Eq. (|76[) at zero magnetic field: 

Me) = \ J2 5(e-e a + ^-E c - J(5 + 1/4)) 

2 

~4ST2 £ <*(e-e«+M-£c-J(S+l/4)) 

+ ^TT ^ <5(e-e« + M-S c + J(5 + 3/4)). 

(138) 

Using the fact that £nq +s+1 = Es+1/2 — E$ and 
e ^si-s+i = E s-i/2 ~ E s , we find 

r^-o^)^-^ ,i39) 

in accordance with the sketch of Fig. [7] 

We can estimate the effect of level fluctuations on the 
TDOS from the qualitative arguments presented above. 
As we have already demonstrated, at zero temperature 
the height of the maximum in the TDOS is given by 
1/(45") for a given realization of single-particle levels. 
Their statistical fluctuations lead to averaging of this es- 
timated result, i.e., the height of the maximum in the av- 
erage TDOS is given by 1/(45) w [1 + (An 2S ) 2 ] J/(2 J*) 
(see Eq. (|117p ). Taking into account Eq. (|119|) and the 
discussion after Eq. (|120[) , we expect that the height of 
the maximum in the average TDOS is of the order of 
[1 + (l//37r 2 )ln(J^/T)]J/(2J^). Similarly to the zero- 
field average spin susceptibility the statistical fluctua- 
tions of single-particle levels result in logarithmic depen- 
dence with temperature of the height of the maximum in 
the TDOS. 

In the presence of a strong magnetic field, b J, 
and with neglect of level fluctutions, the the typical spin 
5 is of the order of JJajJb. Therefore, the width of 
the maximum in the TDOS can be now estimated as 
J 5 ~ J+b/6. The relative height of the maximum be- 
comes 1/45 ~ J5/J+b. Both estimates are in agreement 
with the results of Sec. IVIII Bl 

As was demonstrated in Sec. IVI CI a magnetic field, 
b ^> J, strongly suppresses the effect of level fluctuations 



V(r(e)/Vo 



1 



2 









| -JI2J, 










e, e 2 e 3 



FIG. 7: Sketch of dependence of the tunneling density of 
states on energy at zero temperature. Shaded areas are equal, 
(see text) 

on the spin susceptibility. The same holds for the TDOS. 
The level fluctuations result in relative corrections to the 
TDOS (|T37j) of the order of (S/b) ln(J*&/J<5) < 1. 



X. CONCLUSIONS 

We have addressed here the interplay of charging and 
spin-exchange interactions of electrons in a quantum dot. 
Even within the simple UH framework, this problem 
becomes non-trivial due to the underlying non-Abelian 
action that necessarily requires the evaluation of time- 
ordered integrals. To overcome this obstacle we have 
employed here a technique based on the WNK trans- 
formation. It allows us to obtain exact analytic resuts 
which describe the partition function (Eq. (|57| ) and the 
tunneling density of states (Eq. (|76[0 for an arbitrary 
single-particle spectrum, temperature, Zeeman splitting, 
charging and exchange energies. Our solution (|57[) for the 
partition function reproduces the result obtain ed p revi- 
ously by means of another approach in Refs. p5ll26l]. We 
believe that the approach employed in this paper is more 
manageable for analytic calculations and extensions. Our 
result (|7d| for the TDOS is a generalization of the result 
derived in Ref. [13j to the case of finite spin-exchange 
interaction and Zeeman splitting. 

In the mesoscopic Stoner regime, near the Stoner in- 
stability, S — J <C <5, we have analyzed our general re- 
sults ([57)1 and (fTo) in detail. In particular, we have found 
that in a wide temperature range 6 <C T <C 5 J/ (5 — J), 
the average zero-field spin susceptibility behaves accord- 
ing to the Curie law with a large effective spin which 
depends on temperature logarithmically (see Eq. (|107[) h 
This dependence results from statistical fluctuations of 
single-particle levels in QDs. We have demo nstrated that 
a tiny magnetic field B ~ y / JT(l — J/6)/gfiB is suffi- 
cient to suppress temperature dependence of the aver- 
age spin susceptibility, as well as to diminish a role of 
statistical fluctuations. We have found that enhanced 
spin correlations, resulting in a large total spin in the 
ground state of a QD (in the mesoscopic Stoner regime, 
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S — J <C S) , become apparent as additional (to Coulomb 
blockade) non-monotonic behavior in the TDOS (at high 
temperatures, 5 <C T <C S J / (5— J)). Similarly to the case 
of the spin susceptibility, we have found that magnetic 
field suppresses the spin-related non-monotonic behavior 
of the TDOS. 

To test our results one needs to explore QDs made 
of materials close to the thermodynamic Stoner instabil- 
ity. The long list of such materials includes Co impu- 
rities in a Pd or Pt host, Fe or Mn dissolved in vari- 
ous transition-metal alloys, Ni impurities in a Pd host, 
and Co in Fe grains, as well as new, nearly ferromag- 
netic rare-earth materials J£r— For the closest material 
to the Stoner instability, we are aware of, YFc2Zn2o, the 
parameter J/(S — J) is approximately equal to 16, lead- 
ing to the spin in the ground state of the order of 10. 
Our results for the spin susceptibility can be checked 
by measuring of the total magnetization and the elec- 
tronic spin resonance signal. Although, as is well-known 
(see e.g. Ref. do[), a single Fe impurity in a nearly fer- 
romagnetic material typically acquires an effective spin, 
•Simp, of the order of 10, the temperature behavior of 
the impurity's contribution to the spin susceptibility is 
different: x imp ~ S? mp /(T]nT/T ), where the temper- 
ature To is set by the interactions between the electron 
spins and the impurity^i We speculate that the intrigu- 
ing magnetic behavior observed recently in Pd nanopar- 
ticles capped with different protective systems^ can be 
related to the physics of mesoscopic Stoner regime. Our 
predictions for the TDOS can be manisfested in the non- 
linear current-voltage characteristics of QDs, measured 
in the sequential tunneling regime at temperatures lower 
than the bias voltage. The expected spin-related non- 
mo notonicities may exceed 7—15% and should be sen- 
sitive to the applied magnetic field. 

The approach used here for analytic computation of 
the spin susceptibility and the TDOS in the UH with 
isotropic spin-exchange interaction can be applied for 
other problems. For example, it can be used for the 
study of transport through a QD described by the UH in 
the co-tunneling regime that requires calculation of two- 
particle Green's functions of electrons in QD. One more 
example is provided by the UH with anisotropic (XXZ) 
spin-exchange interaction^ This model is of interest 
since it describes a crossover from the case of isotropic 
spin-exchange interaction with mesoscopic Stoner insta- 
bility and trivial dynamical spin susceptibility to the case 
of Ising spin-exchange interaction without mesoscopic 
Stoner regime but with interesting dynamical spin sus- 
ceptibility Also our approach can be fruitful for investi- 
gation of the same crossover from isotropic to Ising spin- 
exchange interaction realized in QDs with spin-orbit in- 
teraction ^ 33 ' 34 Yet another example is given by the com- 
petition between superconductivity and ferromagnetism 
in the UH with an attractive interaction in the Cooper 
channel.— The method used in the present paper can be 
generalized to study this interplay as well. 
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Appendix A: Wei-Norman-Kolokolov transformation 

In this appendix we demonstrate relation of the non- 
linear transformation (|24[) with the general method by 
Wei and Norman^ and calculate the Jacobian (|26[) of 
Wei-Norman-Kolokolov transformation.— 

The equation for time evolution operator for the time 
dependent Hamiltonian H(t) 



(Al) 



with initial condition U(0) = 1 can be represented by 
a finite product of k exponential operators^ The index 
k is a dimension of the Lie algebra generated by H{t) 
where the Hamiltonian is assumed to be linearly depen- 
dent on the group generators. The general formalism of 
time evolution operator construction is known as Wei- 
Norman method^ 

Consider a time-evolution operator for the system 
with linear realization of dynamical SU (2) symmetry de- 
scribed by the Hamiltonian 



H(t) = -0(t)s 



(A2) 



Then, matrix U(t) is given by time-ordering expo- 
nent p3]). The solution of (TAT]) can be parameterized 
as 

U(t) =cxp(/_(£) S -)cxp(/ 2 (t).s z )cxp(/ + (t).s+) (A3) 

where s ± = s x ±is y and s z are three generators of SU (2) 
group. Functions f-(t), f z (t) and f+(t) satisfy the sys- 
tem of differential equations 



if- = -¥ + + o z f- + \e-f_ 

if, = -0 z -0-f- 
if+ = 



(A4) 



■\9-e-t- 
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with initial condition /_(0) = f z (0) = /+(0) = 0. This 
system of equations can be easily obtained with the help 
of Hausdorf formula. — The solution of the system (|A4|) 
depends on the solution of the single Riccati equation 
(the first equation in (|A4[1 ). Parametrization of (|A4[) by 
three new functions /c* (t) , p(t) defined as 



K+(t) =/_(<), 

K-(t) = -if+(t)e?'®, 
P (t) = -if z {t). 



(A5) 
(A6) 



leads to the Kolokolov representation of the time- 
evolution operator^ 

U(t) = e K+ ^ s ~ e is ° So p(*i) d *i 

x exp (is + J k- (h)e- 1 ti 1 p ^ dt2 dt^j . (A7) 

The initial condition /_ (0) = is translated to n + (0) = 
0. The initial conditions / z (0) = /+(0) = are satisfied 
by construction of the functions k~ (t) and p(t) . There- 
fore, the variables n~ (t) and p(t) are not constrained by 
the initial conditions. 

In order to find the Jacobian of the Wei-Norman- 
Kolokolov transformation, we consider its discrete ver- 
sion. Let us split the interval (0, t) into N — > oo parts 
of length A — > such that t = NA. Also we introduce 
9j = 6{tj), = K ± (^j) and pj = p(tj), where tj = jA 
with j = 0, . . . , N. Then, the solution of Eqs (|AT ]) -(|A"2 j) 
can be written as U(t) = Y\f =0 exp(iA0js). Together 
with Eq. (|A7|) it allows us to find the following map corre- 
sponding to Wei-Norman-Kolokolov transformation (see 
Eq. 05) with p = +): 



K + + K + 
Pj J o J ^ -4-1*74, (A8) 



Oj = Pj -Kj{4 + 4-J- 

The map (|A8|) is supplemented by the initial condition 
Kq = 0. The Jacobian of the map (|A8|) is given as 



(A9) 




3=1 



Finally, taking the limit N —> oo and A — s- 0, we obtain^! 



A)- N exp[ l - j dt'p{t' 



(A10) 



in agreement with Eq. (f2l>)) . 



Appendix B: Integration over k^ v in Eq. (|34[1 

In this appendix we demonstrate how integration in 
Eq. over fields n^ p can be performed. We need to 



evaluate the following functional integral: 



z-v 




ipk p (t p )- J dtk p ^ p (t)e o 
o 



ip J dt' p- p (t') 



if f dtp p (t) 



(Bl) 



Let us start from computation of the 2-point correla- 
tion function 



(<(g«; p (t))o = 



J V[ K ±P] K r(t p ) K -r{t)e- ^ ^ 



(B2) 

In order to evaluate the functional integral, we split the 
interval (0, t p ) into N parts of length A — > such that 
t p = NA. We define K p P j = Kp P (tj), where tj = jA with 
j = 0, . . . , N. Then, the 2-point correlation function (|B2[) 
can be written as (t = nA): 

, P K - P) _ (njLi j dt^^K-'e-w- 1 ^ 

\ K p,N K p,n/0 — 



(n;= 1 /<^%,?])e-^ r - 1 ^ 



Here nf p = (n p ^, Kf p N ) and 



±p 



r = j 



/ill 

1 1 



(B3) 



(B4) 



V o o o 



We remind that k p satisfies the initial condition k p = 0. 



Hence, we find 



(K p -(t p )K- p (t)) = J6(t p -t), 



(B5) 



where the Heaviside step function is defined as 8(0) = 1. 
For obvious reasons, the other 2-point correlation func- 
tions vanish. Next, we can find the following set of iden- 
tities for an arbitrary function f(t) and non- negative in- 
teger k: 

((n p p (t p ) f P dtf(t) K ~ p (t)) ) =k\(jj\tf(t) 



ii 



fdKpdx'e-M l J {x p f*dtf{t)** p ) 



J dztpdzt* e 



\* P \ 2 /J 



(B6) 



We see that the functional integral (|B1|) can be substi- 
tuted by the usual one with complex conjugated variables 
>c p and K p corresponding to k p and hl~ p , respectively. 
Therefore, Eq. (|B1|) becomes 



f / dtp p (t) 



TT f d>t P dK *P c -\x P f/J 
j>=± J n 




ipe""- p k p - tC_ v I 
o 



ip f dt'p- p (t') 



(B7) 
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Integrating over variables x p and x£, we find that inte- 
gral (|B7|) equals 



\p=± J Vp=± 



x / " dte- ipf o dt ' p ^ t '^ 



¥(*+-*-) 



(B8) 



This result yields Eq. ([55]). 



Here, the function £2 (ft) is given as 
L 2 (ft) = 16T 2 sinh 2 - J -^Ru/o 



1 coth § 
cosh ft — 1 



f 



| coth § — u; coth ui 
cosh ft — cosh 2ui 



(C8) 



Using the conditions C(fti,ft.2) = C(ft.2,fti) and 
C(ft, 0) = 0, wc obtain Eq. (p?Tj) in which the function 
L(ft) is related with L2(h) as 



Appendix C: Correlation function C(h\, I12) 

In this appendix we present derivation of the re- 
sults (|92p for the correlation function (pTT|) . Also, we 
derive the result (|101[) . 

The single-particle density of states v (E) has non- 
Gaussian statistics. However, the function V(h) is a 
Gaussian random variable 2 ^ since it involves a large num- 
ber of single-particle levels: max{|ft|, T}/S ^> 1. We re- 
mind that the 2-point correlation function of the single- 
particle density of states is given as 2 ^ 



where 



f 



-R 



u/o 



Ru(x) 
Ro{x) 



sin 2 x 



■ 2 
Sill X 



d sin a 

dx x 



sini 



dt. 



(CI) 

(C2) 
(C3) 



L 2 (h) = L"(h). 



(C9) 



To estimate the function L2(h) at T > we can use 
the asymptotic expression of the function Ru/o & t large 
values of its argument: 



Ru/o{x) = 



j3x 2 



x > 1. 



(CIO) 



Here we recall that (3=1 for the orthogonal ensemble 
and (3 = 2 for the unitary ensemble. Then, at T S 
Eq. (|C8|) becomes 



L 2 (h) 



4sinh 2 (/i/2) 
/3^2 



00 duj_ 

^2 



I coth I 
cosh ft — 1 



1 



§ coth | — wcothw 
cosh ft — cosh 2ui 



(Cll) 



At |ft| <C 1 we expand the right hand side of Eq. (|C11[) 
and find 



Using Eqs (J86J), d90J) and ([CI]), wc obtain 

/ ttTlu\ 



C(h 1 ,h 2 )=T 2 J^R 



In 



In 



sinh 2 (^) 



-In 



x In 



cosh 2 (f 
sinh 2 (^L) 

sinh 2 (^) 



cosh 2 ( £^ZH) 



sinh 2 (^) 
cosh 2 (f) 



where we used the identity 

R(j/o(x)dx = it. 

Let us introduce the function 

d"+ m C(fti,ft 2 ) 



C nm (fti, ft 2 ) 



d n h 1 d m h 2 



(C4) 



(C5) 



(C6) 



with integers n,m ^ 1. Then one can check that the 
following exact relation holds 

Cu(hi, ft 2 ) = i 2 (fti + ft 2 ) - L 2 (hi - ft 2 ). (C7) 



i 2 (ft) 



Cl ft 2 

2/37T 2 ' 



where the numerical constant 



ci 



cfw J 1 w coth w — 1 

^[3 " 



sinh 2 uj 



0.37. 



(C12) 



(C13) 



Hence, using Eq. (|C9[) . we obtain the asymptotic expres- 
sion (EH). 

At |ft| » 1, we rewrite Eq. (jClip as 



i 2 (ft) 



2 sinh 2 (ft/2) J 



/3tt 2 
1 







lftl-2 



I ft. I — 2uj coth a; 



I' 1 !/ 2 du 

UJ 2 

00 dw 

\h\/2 W 2 

2 



2 sinh 2 (ft/2) 
|ft|-2 



cosh(ft) 
|ft| — 2ui coth w 



2sinli (ft/2) 
lftl-2 



,\h\ 



\h\ - 2oj 



^ I 1 " 2 



2 sinh 2 (ft/2) 
ftl 



-.,2u; 



f'2 



(C14) 
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where 



c 2 = — / — 1 - w coth w + / 

Jo ^ Jl 



duj lnw 
sinh 2 uj 



0.43. 



(C15) 

Using Eq. (|C9|) . we obtain the asymptotic expression (|9"2"]l 



from Eq. (|C14|) . 

Equation (|99p implies that mean squared fluctuations 
of the level spacing A can be written as 



(A - S) 2 , <S 2 

^^ = C 22 (0,0)^. 



(C16) 



As follows from Eqs ([C7} and (fCi~2|) . C 22 ( 0,0) 
2c 1 /(3tt 2 . Hence, we obtain Eq. (fTOTj) from Eq. (|C16[) 



The last term in the right hand side of Eq. (|D3|) can be 
described as the contribution from independent Gaussian 
variable v(u) = v(—u) with the correlation function 

V 



v{uj)v{u k ) = -^^{uj ~ u k ) 2 \n( Uj - u k f. (D5) 



Then, we find 



y 3/2 e y/4 



du-j 



v 1 + ^J 



x exp 



n 2 y 2 In 2 




Je"^). (D6) 



Appendix D: Evaluation of x(T,0) in the region II a 

In this appendix we perform evaluation of the averaged 
zero-field spin susceptibility in the leading l ogarithm ic 
approximation: in each order of expansion of x(T, 0) in 
powers of 1/((3tt 2 ) we take into account the term with 
the highest power of In 2 J*/T. 

Let us define Sio(a;,y) = d'E,(x,y)/dx, then at S <C 
T < J* (y > 1) 



x(r,o) 



(J*/J) 2 d 
3T dy 



lnSio(0,y). 



(Dl) 



In order to evaluate lnSio(0,y) we use the replica trick. 
For non-negative integer values of n we obtain that 



[H 10 (0, y )] n = I 



2,3/2^/4 



2 1 vv' 



x exp 



1 n t 



(D2) 



Provided the dominant contribution to the integral (|D2[) 
comes from regions with \uj\ <C -^y, we can expand the 
2-point correlation function (|91[) in powers of Uj and u k 
to the second order whereever it possible. We thus find 



y 2 ln2 1 



+ 2bujU k H — (uj + it fc ) 



Here we should still perform averaging over v(u) with 
the help of Eq. (|D5[) . Introducing Gaussian variable z 

to decouple term fx^=i u j^j m the right hand side of 
Eq. (|D6[) , we rewrite it as 



[z w (o,y)V = y 3n/2 e ny/i 



exp 



nVln2 

2/37T 2 



><n 



^ e -a« J 2 + (c+2 ZV / 6)« J ^ + 7J e v 



( u j). 



3 = 1 



(D7) 



We note that the typical values of Uj contributing to the 
integral in Eq. (|D7[) is of the order of (c + 2zVb)/a. Next 
we introduce the variables Xj = Uj -(c+2zVb)/2a. Pro- 
vided typical values of |c + 2z\fb\j2a <C y/y, the limits of 
integration over Xj are the same as for Uj. Taking into 
account that the correlation function (|D5[) is translation- 
ally invariant, we thus find 



— n n n/2 y 3n/2 e n(y+cy(a-nb))/i „z y 2 

[aio(0,i/)] = 2 "a(«- 1 )/ 2 (a-n6) 1 /2 



-e 2 f 



hi 2 
13 — 



f dz _ z 2 / 2z\Jb(a — nb) + C\/a\ n 
J V™ v (a - nb)^/ay J 



(D8) 



Here, we used the fact that u(a;) is an even function. In 
the limit n — > 0, we obtain 



(uj - Mfe) 2 ln(uj - u k ) 2 , (D3) 



' 2/3^ J 

where we introduce the following parameters for conve- 
nience: 



a = 1 



ny_f 
2/3n 2 V 



In 



y 



2(3tt 2 
2ny 3 / 2 
(3tt 2 



lny H 
In 2. 



16 
3 



(D4) 



hiHi O (0,y) = - 



1 + 



In* 



In 



l + zV21ny/(/37r 2 



where 



X 



du 



-u 2 +v(u) 



(D9) 



(D10) 



As we discussed above, the integral over z in Eq. (|D9I) is 
constrained by the condition | z | \Jh\ y/((3ir 2 ) <C 1 . In this 
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case, the typical values of z contributing to the integral 
are of the order of unity. Our assumption is thus self- 
consistent for y/hxy / (/3tt 2 ) <C 1. Therefore, the integral 
over z in Eq. (|D9[) is proportional to \ny/(3ir 2 <C 1, i.e., 
the integral is small compared to the first term. 

Evaluating In A" to the lowest order in y/{f3iT 2 ) with 
the help of Eq. JDSJ) we find 



Here we introduce the function 



F(*,V)=I rf te-^ ^,^? C ° S( r /2 L . (E2) 



sinh (irt) + cos 2 (7rx/2) ' 



In A" =a\ 



y 



(3tt 2 '"' 
oi = (ln2 + 7-2)/4: 



-0.18. 



(Dll) 
(D12) 



Subtituting this result into Eq. (|D9[) , we reproduce 
Eq. pU7j) with the help of Eq. (|DlT) . To the second order 
in y/(/37r 2 ) we find 



where 



"2 



3 /, „ 7\ 2 3tt 2 19 



(D13) 



In 2 + 7 - 



32 24 



dudv 
-0.039. 



-2( U 2 + ™ +t , 2 )/3 u 2 ln(w 2 )u2 ln(u 2 



(D14) 



Equation (|D13|) demonstrates that the result (|107|) (ob- 
tained by expansion in powers of the correlation func- 
tion <j90j) ) is valid provided the condition y/((3-K 2 ) <C 1 
holds. 



Appendix E: Asymptotic expression of the function 

T(x,y) at y > 1 



which obeys 



?{x + 2k,y) = (-l)\ k \?(x,y), F{2k + l,y) = 

(E3) 

for integer k. At y 3> 1 we find 



~. . sgn[cos(7rx/2)l / ycos 2 (7rx/2) 
V) = -, /„ ; ex P 



(E4) 



y 1 /4^ C os( 7 ra;/2) ™" r V 
y cos 2 (772/2) ' 



x W_i i 

4 ' 4 



where W\^(z) denotes the Whittaker function. However, 
the Whittaker ^-1/4,1/4(2) function is related with the 
error function: 



W_ i i (z) = V^z 1/4 e z/2 [l - erf 



(E5) 



In this appendix we outline derivation of the asymp- 
totic expression (|135[) of the function T{x,y) at y ^s> 1. 
Provided (2n + 1) < |x| < (2n + 3) with integer n ^ 0, 
the function ^(x, y) can be rewritten as follows: 



F(x,y) = -e- v/4 e yx/2 e- yx ' /4 lT(x,y) 



2^(-l) m exp 



m=0 



(N-(2m+l)) a 



Therefore, for y> 1 wc obtain 



F{x,y) = sgn[cos(7rx/2)] exp 



y cos 2 (7rx/2) 



1 — erf 



yy|cos(7ra/2)| 



( E1 ) Using Eqs dETJ and (|E6]) . we find Eq. ([T55]> . 
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